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ABSTRACT 

We have observed a 3° x 3° area centered on the M81/M82 group of galaxies using the Robert 
C. Byrd Green Bank Telescope (GBT) in a search for analogs to the High Velocity Clouds (HVCs) 
of neutral hydrogen found around our galaxy. The velocity range from -605 to -85 km s" 1 and 25 
to 1970 km s _1 was searched for H I clouds. Over the inner 2° x 2° the la detection threshold 
was 9.6 x 10 5 M Q . We detect 5 previously unknown H I clouds associated with the group, as well 
as numerous associated filamentary H I structures, all lying in the range —105 < Vheiio < +280 
km s _1 . From the small angular distance of the clouds to group members, and the small velocity 
difference between group members and clouds, we conclude that the clouds are most likely relics 
of ongoing interactions between galaxies in the group. 



Subject headings: galaxies: clusters: individual(M81) 
galaxies 



1. Introduction 



H I clouds are found around some spiral galax- 
ies at velocities not simply related to the rotation 
of the galactic disk. In the Milky Way there are 
the high velocity clouds (HVCs) both ionized and 



galaxies: interactions — surveys — radio lines: 



neutral found in large complexes and as smaller 
isolated clouds which together cover nearly 40% of 



the sky (Wakker & van Woerden 1997 Lockman 
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et al. 20021. M31 has a similar system (Thilker 
et al. 2004 1 . These clouds most likely result from a 
variety of phenomena. Interactions between galax- 
ies can produce filamentary structures of gas and 
tidal streams like the Magellanic Stream (|Putman" 



et al. 2003 ) and the tidal tails of the M81 group 
flYun et al.[| 19^|Walter eraI1[2002| . Some HVCs 



may be the remnants of galaxy formation that 



are currently being accreted (Mailer & Bullock 



2004), or clumps and filaments of H I associated 



with faint dwarf galaxies ( 


Thilker et al.|2004 


Ben 


Bekhti et al. 


2006 


Begum et al. 


2006 


) . Some 
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clouds may also arise from a "galactic fountain" , 
wherein hot gas created by supernovae expands 
thermally cools, and condenses and falls back into 



the galaxy (see Wakker & van Woerden (19971 



and references therein). Studies of cloud prop- 
erties such as mass, morphology, and kinematics 
are necessary to distinguish between these possi- 
bilities. 

Measurement of many physical quantities is de- 
pendent on an accurate value for distance. The ex- 
act distance of the clouds in our galaxy is very dif- 
ficult to measure, so determination of their proper- 
ties has been difficult. Many attempts have been 
made to measure the distances of clouds in the 



Milky Way (see, e.g. 


Schwarz et al. 


(199f 


) ) ; |van 


Woerden et al. ( 


1999 


); 


Ben Bekhti et al. 




>006)). 



However, these efforts have succeeded only in plac- 
ing limits on cloud distances. 

Other galaxies may also have these high veloc- 
ity clouds. For extragalactic sources, distances 
from a H I cloud to its associated galaxy can 
be determined to a high degree of accuracy, so 
their properties may be studied more readily. We 
observed the M81/M82 cluster with the GBT in 
order to map its neutral hydrogen cloud popula- 
tion. The group contains optically bright galaxies 
M81, M82, NGC 3077, and NGC 2976, which show 



the remnants of strong interactions (Walter et al. 



2002), as well as over 40 dwarf galaxies. See Ta- 
ble 1 for a summary of galaxies within our field of 
observation. The M81 group has been the subject 



of many studies including, e.g., Roberts ( 


1972), 


Appleton et al. (1981), Lo & Sargent 


(1979 


1, Yun 


et al. 


(1993 


I, and Boyce et al. (2001 


1 . Lockman 


et al. 


(2002 


) also found diffuse neutral hydrogen 



emission in the direction of the M81 group. They 
find locations where the emission is likely due to 
the interactions of M81/M82 but also find lower 
level emission away from known galaxies. 

In addition to our velocity range search of -605 
to -85 km s _1 and 25 to 1970 km s , we have 
achieved a lower mass detection threshold than 
previous studies. Our observations have a mass 



sensitivity 50-100% lower than the studies by Ap- 



pleton et al. (1981) and Boyce et al. (2001 ), due to 



our high velocity resolution and the excellent sen- 
sitivity of the GBT. Our observations cover a large 
angular area in order to make a complete study 
of the H I cloud population and properties. The 
beam size of the GBT at L band is well matched 



to the expected angular size of these objects, so 
our study provides new insight into their struc- 
ture. The study of IBoyce et al.l (120011) using HI- 



JASS data covered a larger angular area (8° x 8°) 
than our observations, but with a much smaller 
velocity range search (-500 to 500 km s _1 ) and a 
coarser (26 km s"" 1 ) velocity resolution. 

In addition to distinct clumps of hydrogen, 
galaxies are found to have filamentary structures 
of H I gas and tidal streams associated with in- 



teractions ( 


Yun et al. 


1993 


Walter et al. 


2002). 


Begum et al. 


( 


2006 


) also found clumps and fila- 



mcnts of H I associated with faint dwarf galaxies. 
The M81 group is notoriously rich in large-scale H 
I filaments, a clear signature of interaction between 
the galaxies. We include in our study a search for 
such filamentary structures and smaller scale fila- 
ments that might provide additional information 
about interactions. 

H I clouds are also related to dwarf galaxies. 
In a study of nine Local Group dwarf galaxies, 



Bouchard et al. ( 2006 1 found that the concentra- 
tion of HVCs was markedly enhanced near dwarfs, 
and that there were H I clouds near (in projection) 
most of the dSphs and dlrr/ dSphs of the LG, al- 
though they were typically offset from the optical 
center. |Blitz & Robishaw ( 2000[ ) found that many 



H I envelopes of dwarf galaxies in the local group 
were very similar to extragalactic H I clouds in 
mass and size, and although most neutral hydro- 



gen clouds are thought to contain no stars ( Siegel 
et al. 20051, they found that two previously cat- 



aloged clouds actually harbored dwarf galaxies. 
These results suggest a strong relationship be- 
tween H I clouds and dwarf galaxies, although 
the nature of the connection has yet to be fully 
elucidated. Perhaps some H I clouds enshroud 
low surface-brightness galaxies, or are ejected from 
dwarf galaxies through other mechanisms. 



Bouchard et al. ( 2006 1 concluded that ram pres- 



sure or tidal stripping are not likely to be the cause 
of the observed offset of H I clouds from their as- 
sociated galaxies. An alternative explanation is 
that the gas could be blown out by star forma- 
tion. Recently, Walter et al. ( 2002 1 have shown 



that episodes of star formation in the dwarf galaxy 
NGC 3077 have disrupted and ejected the gas out- 
side of the region where the stars are optically 
identified. Our observations cover a large angu- 
lar area around this galaxy (as well as others in 
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the group) to search for these blown out compo- 
nents. In addition, these larger angular regions 
enable us to distinguish foreground components in 
our galaxy and smaller angular size components of 
dwarf galaxies in the group. 

In sections |2] and |3] we summarize the observa- 
tions and data reduction. In section 2] we describe 
methods used to interpret our observations. In 
section [5] mass calculations for known objects in 
the group are presented. In section|6]we detail new 
detections, including H I clouds and filaments, and 
their properties. In section [7J our new detections 
are compared with a numerical simulation of the 
M81/M82 group. Finally, in section [8] we discuss 
the implications of our results and future work. 

2. Observations 

We observed the M81/M82 group with the 
100m Robert C. Byrd Green Bank Telescope 
(GBT) of the NRACQon 12-16 June, 2003. The 
full group extends over a 40 x 20 degree area (right 
ascension range 7 to 11 hours, declination range 
58° < S < 78°), but the brightest galaxies (M81, 
M82 and NGC3077) are all within a 2° x 2° region 
centered on M81, which was where we concen- 
trated our observations. The 2° x 2° region was 
observed by moving the telescope in declination 
and sampling every 4' at an integration time of 
5 seconds per sample. Strips of constant declina- 
tion were spaced by 2/7. A larger field was then 
sampled at 4' spacing and an integration time 
of one second per spectrum. The data discussed 
here cover an area 3° x 3° centered on M81 to 
include NGC 2976 and the extended H I emission 
associated the group. 

The typical system temperature for each chan- 
nel of the dual-polarization receiver was 20 K. The 
GBT spectrometer was used to cover a velocity 
range -3172 < Vheiio < 7376 km s _1 . Spectra were 
smoothed to an effective channel spacing of 1.29 
km s . Over the inner 4 square degree region, 
the resultant typical RMS noise is 18 mK; outside 
the central region the noise in a channel rises to 
ss 51 mK. The instrumental parameters are sum- 
marized in Table 2, and Table 3 gives a summary 
of observations. 



3. Data Reduction 

The archived GBT data were reduced in the 
standard manner using the GBTIDL^] data reduc- 
tion package. 

In order to match our velocity resolution to 
the expected linewidths of H I clouds in the 
group, we smoothed the data to contain 2048 
channels with a channel spacing of 24.4 kHz, cor- 
responding to a velocity resolution of 5.2 km s _1 . 
A reference spectrum for each set of scans was 
made using the first and last scan of the set. 
The reference spectrum was then used to perform 
a (signal-reference)/reference calibration of each 
integration of each scan. The calibrated spec- 
tra were scaled by the system temperature, cor- 
rected for atmospheric opacity and GBT efficiency. 
We adopted GBT efficiency equation (1) from 
Langston & Turner (20071. We assumed a zenith 
atmospheric opacity To = 0.009. The GBT obser- 
vations were performed without Doppler tracking. 
Within GBTIDL the velocities were shifted to the 
heliocentric frame. 

After amplitude calibration, the spectral band- 
pass variation was removed by subtraction of a me- 
dian filtered baseline with 1300 km s _1 width (6.0 
MHz width). The median filter baseline subtrac- 
tion process is summarized here and is described 



with example figures in Langston & Turner ( 2007 1 



The baseline was subtracted for each polarization 
separately, in order to remove the Intermediate 
Frequency (IF) electronics gain variations that are 
unique to each signal path. The baseline subtrac- 
tion process also removes gain and sky variations 
between the time of observation and the time of 
observation of the reference scan. The baseline 
value is computed for each individual channel of 
the spectrum by computing the median value of all 
channels within a range of channels centered on 
the individual channel. The median width must 
be small enough to preserve the gain variations, 
yet large enough to reject spectral lines within the 
median range. After computation, the baseline is 
subtracted from the original spectrum. 

Calibrated GBTIDL 'keep' files were processed 
using idlToSdf its [^] This program is used to 



lr The National Radio Astronomy Observatory (NRAO) is a 
facility of the National Science Foundation operated under 
cooperative agreement by Associated Universities, Inc. 



2 Developed by the National Radio Astronomy Observatory; 

documentation at http://gbtidl.sourceforge.net 
3 Developed by Glen Langston of NRAO; documentation at 

http: / /wiki.gb. nrao.edu/bin/view/Data/IdlToSdfits 



3 



flag data contaminated by radio frequency inter- 
ference (RFI) and convert the data to the form 
required for input into the NRAO AIPS^] pack- 
age. The frequency range observed was relatively 
free of RFI and less than 0.2 % of all spectra were 
adversely affected by RFI. The spectra exhibit- 
ing RFI values were identified by tabulation of the 
RMS noise level in channels free of neutral hydro- 
gen emission. In order to remove RFI from the 
data, idlToSdf its was run on each set of scans, 
and the resulting RMS noise values were exam- 
ined. Spectra that showed high values of RMS 
noise across many channels were flagged and re- 
moved. After calibration, the observations were 
imported into the AIPS package. Observations 
were gridded using the AIPS task SDIMG, which 
also averages polarizations. 

The beam size was found by fitting a two di- 
mensional Gaussian to the continuum emission 
from M82. The fitted beam was 10.1'± 0.6' x 
9.4'± 0.5', at a position angle of 54° ± 32°. In 
order to convert to Jy/beam, we observed the cal- 
ibration source 3C286. Ott et al. ( 1994 1 provide 



In their study of H I clouds in M31, Thilkcr 



et al 
7=r 



a model for the flux density of 3C286 relative to 
other radio sources. For 3C286 at 1.41 GHz the 
flux density is 14.61 ± 0.94 Jy. Based upon the 
ratio of 3C286 to M82 intensities, the flux density 
of M82 is 7.17 ± 0.46 Jy. Our value for M82 is 
consistent with 6.2 ± 0.2 Jy from Condon et al. 
( 1998]) . We determined the calibration from from 
K to Jy by producing a continuum image of M82 
in the same manner as the spectral line images. 
The peak M82 brightness temperature (K) was 
measured relative to laboratory measurement cal- 
ibration noise diodes. The calibration factor was 
determined by the ratio of the measured bright- 
ness temperature divided by flux density of M82 
(Jy). Including all corrections for the GBT effi- 
ciency and the mapping process, the scale factor 
is 0.47 Jy/K. 

4. Data Interpretation 

For all estimates of H I masses, we adopt 
the M81 distance, D = 3.63 ± 0.34 Mpc from 
Karachentsev et al. (20041). M81 is at a distance 



between M82 and NGC3077, 3.52 and 3.82 Mpc 
respectively. 



(20041 found clouds mostly within 150 km 
of the systemic velocity of M31. Filaments 
were found within 80 km s" 1 of the galaxy's ve- 
locity. If these values can be interpreted as a 
prediction for H I cloud locations, we expect to 
find clouds with velocities within 150 km s _1 rel- 
ative to the main galaxies in the group. For this 



paper we use galaxy velocities from Yun ( 1999 1 



For completeness, spectral maps were made in a 
wide velocity range, from -605 to 1970 km s _1 . 
This extends the search to more than ± 500 km 
s _1 beyond the systemic velocity of any group 
galaxy. The 7a detection threshold for H I mass 
was 9.6 x 10 5 M Q . Spectral maps were visually in- 
spected for possible new clouds. Objects that were 
spatially distinct from group galaxies and tidal 
streams in a range of velocity were considered to 
be cloud candidates. Positions of candidates were 
measured by fitting a 2-dimensional Gaussian with 
the AIPS tasks JMFIT and IMSTAT to the peak 
emission. For each H I cloud candidate, spectra 
were obtained from ISPEC, and H I column den- 
sity maps for the velocity range in which the can- 
didate was detected were made with MOMNT. 

Our search area of 3° x 3° centered on M81 
corresponds to 194 x 194 kpc 2 . Based on studies 
of the Milky Way and M31, which found clouds 
anywhere within 1-50 kpc of the center of a galaxy 



(Thilker et al. 2004 Wakker & van Woerden 1997 1, 



we expected that H I clouds should be located 
within 50 kpc of associated galaxies. Therefore, 
we should detect all of the H I clouds around M81. 
For M82, NGC 3077, and NGC 2976, our field of 
view excludes some of the regions where clouds 
might be found; as such, our detections for those 
galaxies will place lower limits on the H I cloud 
population. iThilker et al.l d2004| found that fil- 



4 Developed by the National Radio Astronomy Observatory; 
documentation at http://www.aoc.nrao.edu/aips 



amentary structures were limited to the inner 30 
kpc of the galaxy, so we expect to detect all such 
structures associated with the M81 group. 

The boundaries of regions containing clouds 
were determined by inspection of the spectral line 
images. The boundaries of the region averaged 
into the H I intensity profile were chosen to sep- 
arate the individual components from the other 
galaxies. However all of the newly identified clouds 
are in close proximity to the interacting members 
of the galaxies. For each profile, the region was 
set using a box approximating the area of neutral 
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hydrogen most clearly associated with the object 
in question and not involved in interaction. 
The spectral line image is available on-lincQ 

5. Masses of Known Objects 

There is some overlap in velocity space between 
the M81/M82 group and the Milky Way and over 
this velocity range we can not discriminate be- 
tween local and distant emission. To estimate 
the total H I mass of neutral hydrogen for the 
M81/M82 group we replaced the velocity channels 
contaminated by Milky Way and local high veloc- 
ity clouds with interpolated values computed from 
spectral line observations at -85 and +25 km s _1 . 

The column density was calculated for the en- 
tire field observed and each new H I cloud. Mo- 
ment maps showing the column density are pre- 
sented in Figures [T] and [3] For the observed field 
and major galaxies, the flux values at velocities 
from -250 to 340 km s _1 were summed. For new 
objects, the flux in the velocity range over which 
each object detected was summed. The flux vs. 
velocity profile of the object was obtained from 
the MPS task ISPEC. 

5.1. Total Mass over the Field 

Neutral hydrogen was detected over the velocity 
range -250 to 340 km s _1 , and these channels are 
summed and shown as contours in Fig ure [T] The 
integrated column density map from Yun et al. 
( 1994 ) is shown as the grey scale image in Figure [lj 
The locations of each major galaxy, dwarf galaxy, 
and H I cloud detection are indicated. 

The spectrum of the entire field used to com- 
pute the neutral hydrogen mass of the group is 
shown in Figure [2] The computed mass was 10.46 
± 2.86 x 10 9 M Q . For this mass and all following 
mass estimates, the error estimates are la. The 
dominant error contribution is the uncertainty in 
the absolute calibration. Our total mass value is 
greater than the values found by |Appleton et al. 
(1981) and Yun (19991 . Our study is more sen- 



sitive to faint, extended emission over the entire 
field than either of these previous studies. 



5 The spectral line image cube 
|http: //www.nrao.edu/astrores/m81| 



available at 



5.2. Major Galaxies 

The H I spectra for the major galaxies are pre- 
sented in Figure [2] Neutral hydrogen masses of 
known galaxies were calculated to confirm the ac- 
curacy of our mass calculation methods. Our cal- 
culated values for Mjji are in agreement with the 
literature; see Table 4. The dwarf galaxy BK IN 
( Huchtmeier fc Skillman|[l998 1 is also detected by 
our observations in the velocity range ~ 560 to 580 
km s _1 . BK IN has a velocity outside the range 
shown in Figure [2j 

6. New Identifications 

The newly detected clouds were identified by 
individual examination, by eye, of all spectral line 
images for this region in the velocity range -605 to 
-85 km s" 1 and 25 to 1970 km s" 1 . We detect five 
neutral hydrogen clouds associated with the M81 
group of galaxies. See Tables 5 and 6 for a sum- 
mary of detections. Three of the clouds are clearly 
associated with a particular galaxy, and two are 
located between two known galaxies. All of the 
clouds have a velocity close to that of at least one 
of their parent galaxies. Clouds were only detected 
between -105 and 280 km s , and the maximum 
difference between a cloud and associated galaxy 
was 119 km s _1 . 

Figure [3] shows H I column density maps of 
each detection. For these maps, only the velocity 
range in which the cloud of interest was detected 
is summed into the integrated intensity map. One 
channel on each end of the detected range was ex- 
cluded from the integrated intensity map in or- 
der to increase the signal to noise ratio and high- 
light the cloud. Additionally, all of the clouds are 
within 35 kpc (in projection) of the group. How- 
ever, this value should be taken as an estimate, 
as the exact distance to each galaxy in the group, 
and more importantly, the distance to each cloud, 
is uncertain. In any case, the approximate dis- 
tance values are consistent with detected distances 
of clouds near M31 ( |Thilker et al]|2004 l. 

Figure [3] shows the spectra of each H I cloud. 
The mass of each cloud was measured by sum- 
ming the flux densities measure for the angular 
regions over the velocity range show in figure [4] 
The masses of the detected clouds are of order 10* 
Mq. In order to measure the systemic velocity 
and velocity dispersion of each cloud, a Gaussian 
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function was fit to the average intensity profile. 

Figure [5] shows position- velocity (PV) diagrams 
for each new H I cloud. Whereas clouds may 
be difficult to distinguish from the extended H 
I emission associated with major galaxies and tidal 
streams in the column density maps, the H I spec- 
tra and position-velocity diagrams show the dis- 
tinct nature of these new clouds. 

We have compared the GBT observations with 



the VLA images from Yun et al. ( 1993 ) . Below, we 



describe the H I distribution found in the VLA im- 
ages at the locations of each cloud. The combined 
GBT and VLA images are presented in |Langstoh] 
et al., (in prep I, along with a description of the 



process of combining GBT and VLA images. 

6.1. Cloud 1 

Cloud 1 is located to the north-west of M82, 
separated by approximately 31 kpc in projection. 
It is clearly separated from M82. It is detected in 
the velocity range 132 to 204 km s -1 , with peak 
emission at 165 km s _1 . We calculate a mass of 
1.47 ± 0.35 x 10 7 M Q for Cloud 1. 

Cloud 1 lies outside the region imaged by Yun 



et al. (19931, and no VLA images are currently 
available for this area. VLA observations are 
scheduled to supplement the previous work. 

6.2. Cloud 2 

Cloud 2 sits immediately to the south-east of 
NGC 3077, at a distance of 23 kpc. It is a distinct 
clump of H I situated within a filamentary struc- 
ture that extends from NGC3077 south and west, 
encompassing the dwarf galaxies BK5N, IKN, and 
K64. The cloud is in the upper range of the fil- 
ament velocity. The morphology of the filament 
to the west of the cloud suggests that it may be 
plunging into NGC3077. Alternatively, the cloud 
may be gas blown out of the galaxy by star for- 
mation. 

Cloud 2 is detected in the range -162 to -69 km 



with peak emission at -105 km s 



We cal- 



culate a mass of 2.25 ± 0.49 x 10 7 M m for Cloud 



Cloud 2 is very weakly visible in the VLA im- 



ages of Yun et al. ( 1993 1. Clearly the VLA images 



are missing diffuse emission associated with this 
cloud. Other regions have similar strength in the 
VLA images, but are significantly weaker in the 



GBT images. A weak secondary peak in the GBT 
image, at 10 /l 02 m 12.3 s , 68°20'34", is also visible 
in the VLA image. 

This cloud and the associated filament may be 
associated with the clouds in the same vicinity de- 



tected by Brinks et al. (2007) using the VLA. We 



will seek to confirm the detections by |Brinks et al.| 
( 2007 ) with our scheduled VLA observations of the 



region. 

6.3. Cloud 3 

Cloud 3 is located 27 kpc to the north-east 
of NGC 2976. It is distinct from NGC 2976. 
Cloud 3 overlaps with the velocity range of the 
Galaxy, making properties difficult to determine 
accurately. Cloud 3 is detected in the range -17 
to 86 km s _1 , with peak emission at 11 km s _1 . 
We calculate a mass of 2.67 ± 0.65 x 10 7 M for 
Cloud 3. 

Cloud 3 is clearly resolved in the GBT image 
and only barely detected above the background 



noise in the VLA images of Yun et al. ( 1993 1. 



6.4. Cloud 4 

Cloud 4 is located to the west of the H I bridge 
between M81 and M82. It is clearly bound to the 
H I bridge, but is a distinct clump. Cloud 4 also 
overlaps with the velocity range of the Galaxy. 
Cloud 4 is detected in the range -19 to 111 km 
s _1 , with peak emission at 72 km s _1 . We calcu- 
late a mass of 8.37 ± 1.75 x 10 7 M for Cloud 
4. 

Although Cloud 4 is located equidistant from 
M81 and M82 at a projected distance of 33 kpc 
from both galaxies, it is most likely associated 
with M81, as it is much closer to M81 in veloc- 
ity. 



In the H I images of Yun et al. (1993), cloud 



4 has the most compact structure of our newly 
detected clouds. The peak in the VLA emission in 
this region is very close to the center of the GBT 
detection. 

6.5. Cloud 5 

Cloud 5 is to the south of M82, located approx- 
imately halfway between M81 and M82 (separated 
by 18 kpc from both galaxies). Cloud 5 is only dis- 
tinct from M82 in a small range of velocities. It 
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is unclear whether the cloud is actually a separate 
feature from M82; it appears as a red shoulder in 
the H I profile. In order to calculated the mass of 
the cloud a Gaussian was fit to the H I spectrum 
in the peak range where the cloud was detected. 
The cloud is detected from approximately 267 to 
302 km s^ 1 , peaking at approximately 280 km s" 1 . 
Since the H I profile of the cloud is blended with 
M82, the calculated mass is difficult to determine. 
The results of the Gaussian fit give a mass of 0.69 
± 0.27 x 10 7 M G . 

Cloud 5 is a large diffuse region in the maps of 



cold structure (Toomre & Toomre 19721. This 



Yun et al. ( 1993 ), only a few sigma above the noise 



in the VLA images. 

6.6. Other Structures 

All H I cloud detections have some level of H 
I structure surrounding them, connecting them to 
the associated galaxy. In particular, Cloud 2 is a 
clump within a pronounced filament to the south- 
west of NGC 3077. This cloud/filament structure 
may be associated with the dwarf galaxies BK5N, 
IKN, and K64, as they are close in projection. The 
velocity of K64 is known, at (vkg4 — -18 km s _1 ), 
however, the detected velocities of the cloud and 
filament are within about 90 km s _1 of VK64, so 
association is certainly possible. 

Cloud 4 appears to be quite tightly bound to 
the H I bridge connecting M81 and M82. 

7. Comparison With Numerical Simula- 
tion 

A useful insight on whether the newly discov- 
ered H I clouds are part of the system of tidal 
debris produced by the recent collisions involving 
M81, M82, and NGC 3077 can be obtained by 
comparing their observed locations and kinemat- 
ics with the results of a numerical simulation of 
the tidal interactions in this group. Such a com- 



parison can be made with the simulation of Yun 



(1999) 



The simulation by Yun ( 1999 1 is primarily an 



illustrative model in that only the velocity field of 
tidal debris is mapped using a limited three-body 
calculation, rather than employing a full N-body 
calculation. On the other hand, this approach has 
been shown to be effective in tracing the location 
and kinematics of tidally disrupted H I disks as the 
H I emitting clouds originate from a dynamically 



simulation suggests that nearly all of the observed 
H I features in the M81 group are consistent with 
being part of the tidally driven structures result- 
ing from the collisions involving all three galaxies 
during the last 200-300 Myrs - see |Yun| ( fl999"t for 
more detailed discussions. Specifically, the loca- 
tions of the Clouds 1, 2, 4, and 5 are coincident 
with the locations of tidally sprayed material from 
the individual disks. While a chance coincidence 
cannot be ruled out, this comparison suggests that 
the origin of these H I clouds are tidally disrupted 
disk material, rather than a kinematically distinct 
satellite in the group. A more sophisticated, self- 
consistent future numerical simulation will offer an 
improved understanding on the nature of these H 
I clouds. 

8. Conclusions 

We have detected 5 new H I clouds in the 
M81/M82 galaxy group, all with associated fila- 
mentary structures of H I gas. We do not find 
a population of clouds at velocities between -605 
and -105 km s _1 or between 280 and 1970 km 
s _1 . All of the objects have properties similar to 
those of clouds previously found in our galaxy and 
other nearby galaxies. The newly detected clouds 
have small distance offsets from larger neighboring 
galaxies and also have small velocity offsets. 

We have also measured the mass of group galax- 
ies, obtaining values in agreement with previous 
measurements, and calculated masses for our new 
objects. The detected H I clouds have masses 
ranging from 0.69 to 8.37 x 10 7 M . These cloud 
masses are larger than the recently discovered 
clouds near M31, and further high angular reso- 
lution images with greater sensitivity will likely 
yield detection of lower mass clouds. 

The region was searched for clouds from -605 
to -85 km s _1 and 25 to 1970 km s^ 1 . With a 7a 
detection threshold of 9.6 x 10 5 M Q , clouds were 
found only within 120 km s -1 and 35 kpc in pro- 
jection of the galaxy with which they are asso- 
ciated. Therefore, it is likely that these clouds 
are part of the current interactions ongoing in the 
galaxy group. 

Previous studies of H I clouds have been focused 
on the Milky Way and M31. Both of these galax- 
ies are in a fairly relaxed state, and not currently 
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undergoing strong interaction with other galaxies. 
The M81/M82 group, in contrast, is obviously in- 
teracting. Therefore, the study of H I clouds and 
other structures in the group is very important 
in order to understand the full lifetime of these 
structures in relation to galaxy interactions. Our 
results are inconsistent with models of primordial 
H I clouds falling into the cluster. Perhaps this 
conclusion may be extended to the HVCs in the 
Milky Way; however, the clouds must be studied 
further to make this comparison. 

Numerical modeling of ongoing interactions in 
the group are in progress ( |Chynoweth et al.|2008|>. 
Simulations will improve on the model of Yun| 
(1999) by utilizing a fully self-consistent N-body 
model of the galaxy group, with particular em- 
phasis on reproducing the isolated H I structures 
observed. Follow-up VLA observations our newly 
detected clouds are scheduled, and we are prepar- 
ing a description of combined VLA and GBT im- 
ages of these clouds. These higher angular resolu- 
tion interferometer images will be combined with 
GBT data to facilitate numerical modeling. Ob- 
servations at higher spatial resolution will also be 
useful in order to probe the fine structure of these 
objects and detect objects with a smaller angular 
size. 

Facilities: GBT. 
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Tabic 1: Known Galaxies in the M81 Group a 



Name 


a b 


S b 




Type d 


D b 

^ may 


D b 


Inclination 




(J2000) 


(J2000) 


(km s _1 ) 








° E of N 


M 81 


09 ft 55 m 33.5 s 


69°03'06" 


-34 


Sb 


26.9 


14.0 


147 


Ho IX 


09 h 57 m 32.4 s 


69°02'35" 


46 


Irr 


2.5 


2.0 




BK 3N 


09 ft 53 m 48.5 s 


68°58'09" 


-40 


Irr 


0.5 


0.4 




A0952+69 


09 h 57 m 29.0 s 


69°16'20" 


100 


Irr 


1.8 


1.6 




KDG 61 


09 /l 57 m 02.7 s 


68°35'30" 


-116 


Sph 


2.4 


1.4 




M 82 


09 h 55 m 53.9 s 


69°40'57" 


203 


Irr 


11.2 


4.3 


-114 


NGC 3077 


10 h 03 m 21.0 s 


68°44'02" 


14 


Irr 


5.4 


4.5 


-142 


Garland d 


10 h 03 m 42.7 s 


68°41'27" 


50 


Im 








BK lN f 


09 fc 45 m 14.3 s 


69°23'23" 


571 


S 








FM 1 


09 h 45 m 10.0 s 


68°45'54" 




Sph 


0.9 


0.8 




BK 5N 


10^04 m 40.3 s 


68°15'20" 




Sph 


0.8 


0.6 




IKN« 


10 h 08 m 05.9 s 


68°23'57" 




Sph 


2.7 


2.3 




NGC 29768 


09 h 47 m 15.6 s 


67° 54' 49" 


3 


Sm 


5.9 


2.7 


-119 


U5423 g 


10 h 05 m 30.6 s 


70°21'52" 


348 


BCG 


0.9 


0.6 




KDG 648 


10 /l 07 m 01.9 s 


67°49'39" 




Sph 


1.9 


0.9 




KK 77 


oq^so*" io.o s 


67°30'24" 




Sph 


2.4 


1.8 




IC 25748 


10 h 28 m 22.4 s 


68°24'58" 


57 


Sm 


13.2 


5.4 


50 



"Within observed area 
b From 



Karachcntscv ct al. 



(2001) 



Yuii] ([1999} for major galaxies, Karachcntscv ct al. ( 2004 1 for others 



c From 

''From Karachentsev et al 



I 2002 I 



FromjXpplcton et al. ||1981| 



^From |Huchtmeier fc Skillmanj ( |1998 i. 
9 May not be contained entirely within our search area 



Table 2: Parameters of the Robert C. Byrd Green Bank Telescope System 
Telescope: 

Diameter 100 m 

Beamwidth (FWHM) 9.1 ' 

Linear resolution 2.7 Dm pc kpc 

Aperture efficiency ~ 0.69 

Receiver: 

System Temperature ~ 19.5 K 

Spectrometer: 

Bandwidth 50 MHz (10550 km s" 1 ) 

Resolution, Hanning Smoothed 6.1 kHz (1.29 km s _1 ) 
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Table 3: Observations Summary 



Area: 

a range (J2000) 08 ,l 25 m 31.9 s - 11°25'31.9" 

8 range (J2000) 67 /l 33 m 25.7 s - 70°33'25.7" 

Observations: 

Center Frequency (MHz) 1410 

Bandwidth (MHz) 50 

Channel Width (kHz) 24.4 

Velocity Resolution (km s _1 ) 5.2 

Integration time (hours): 

Inner Region 3.88 

Expanded Area 6.75 

Typical RMS noise (mK/channcl): 

Inner Region 18 

Expanded Area 51 

Sensitivity to H I (la): 

N H (x 10 17 ) 

Inner Region 2.5 

Expanded Area 7.0 

M H (x 10 5 Mq) 

Inner Region 1.4 

Expanded Area 3.9 



Table 4: Calculated Galaxy Masses 



Galaxy 


Mass 
(xlO 9 Af ) 


Mass from Yun (1999) 


Mass from Appleton et al. (19811 






Field of Study 

M81 

M82 

NGC3077 
NGC2976 


10.46 ± 2.86 
2.67 ± 0.55 
0.75 ± 0.16 
1.01 ± 0.21 
0.52 ± 0.11 


5.6 

2.81 ± 0.56 
0.80 ± 0.16 
0.69 ± 0.14 


5.4 

2.19 ± 0.22 
0.72 ± 0.07 
1.00 ± 0.10 
0.16 ± 0.02 



Table 5: H I Cloud Detections: Observed Properties 



Cloud 


Association a 


a 


5 


^^assoc 






(J2000) 


(J2000) 


(kpc) 


1 


M82 


09 ft 50 m 07.2 s 


69°55'56" 


31 


2 


NGC3077 


10' l 02 m 50.7 s 


68°19'49" 


23 


3 


NGC2976 


09 fl 52 m 44.9 s 


68°12'50" 


27 


4 


M81 


10 h 01 m 45.0 s 


69°16'31" 


33 


5 


M82 


09 h 55 m 06.3 s 


69°22'05" 


18 



"Most likely association, based on positional and velocity proximity 
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12 3 4 




1010 05 00 09 55 SO 45 40 

RIGHT ASCENSION (J2000) 



Fig. 1. — The contours show the GBT H I integrated column density map of the field observed. Contours 
are overlaid at 1500, 3000, 7500, 15000, 30000, 75000, 150000 and 300000 kJy/beam x km/sec. The velocity 
range included is -250 to 340 km s -1 . The integrated column density map from the VLA observations of 
Yun et al. (19941 is shown in greyscale over the range -400 to 4000 kJy/beam x km/sec. 
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Field of view 






10 








5 

















-5 


i i i L i i 







-400 -200 200 400 

Fig. 2. — H I intensity profiles of major group galaxies, used to calculate H I mass. From top left: M81, 
M82, NGC 3077, and NGC2976. Also included is the H I spectrum of the entire area observed (bottom left). 
Note that the velocity range from -85 to 25 km s _1 is interpolated to remove foreground gas 
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1 04 02 



52 50 



Field of View (-250 to 340 km s" 1 ) Cloud 1 (160 to 185 km s" 1 ) 




Cloud 2 (-115 to -80 km s" 1 ) 




Cloud 4 (75 to 95 kms" 1 ) 




Cloud 3 (5 to 75 km s" 1 ) 



1 

-® 








l7?V~T7 loud 5 - 


JJGC 307 

1 1 


1 1 


1 1 1 1 1 



Cloud 5 (270 to 320 kms" 1 ; 



Fig. 3. — Contour maps of each H I cloud candidate, overlaid with H I column density maps of the velocity 
range in which each cloud is detected. The greyscale is to 12.5 kJy/beam x km/sec. Contour levels are 
65 Jy/beam x km/sec x (10, 20, 50, 75, 100, 200, 500, 10000). The top left-hand panel shows the zoom 
regions for the contour plots; the upper box includes Clouds 1, 4, and 5 and the lower box includes Clouds 
2 and 3. Below each map, the velocity range summed into the integrated intensity map is indicated. 
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Fig. 4. — H I profiles for all cloud detections. The solid line indicates the range that was summed in order 
to calculate the cloud mass. For Cloud 5, the data (dashed line) along with Gaussian fit for Cloud 5 (solid 
line) are plotted. Note that the velocity range from -85 to 25 km s _1 is interpolated to remove foreground 
gas 
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Cloud 1 

— i 1 

+ 69.59.00.1 74 



Cloud 2 




58" 
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h 00 n 

or 
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04" 
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l 1 1 — avi 1 1 1 1 




63.19.59.946 
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Cloud 3 
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-60 -20 20 60 100 140 
Velocity (km/sec) 
Cloud 5 
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Velocity (km/sec) 




09 n 59 
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Velocity (km/sec) 

Fig. 5. — PV diagrams for all H I cloud detections. Each cloud is located approximately in the center of the 
diagram, marked by pointers. Each diagram extends ± 100 km s" 1 from the cloud's peak velocity. The PV 
cut runs through the center of the clouds in declination, and the diagram is centered on the cloud in right 
ascension. The greyscale runs from -0.02 to 0.2 Jy beam -1 . The parent galaxy for each cloud are also visible 
in the PV diagrams. 
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Table 6: H I Cloud Detections: Derived Properties 



Cloud 


Peak Flux Density 


Integrated Flux 




_ a, b 

a v 


Vfl-/-V assoc 


Mass 




(Jy) 


(Jy km s _1 ) 


(km s _1 ) 


(km s _1 ) 


(km s^ 1 ) 


(xlO 7 M ) 


1 


0.11 ± 0.08 


0.91 ± 0.08 


165 


50 


38 


1.47 ± 0.35 


2 


0.10 ± 0.08 


1.39 ± 0.09 


-105 


55 


119 


2.25 ± 0.49 


3 


0.12 ± 0.08 


1.65 ± 0.10 


11 


82 


8 


2.67 ± 0.65 


4 


0.30 ± 0.08 


5.18 ± 0.28 


72 


28 


106 


8.37 ± 1.75 


5 


0.07 ± 0.08 


0.43 ± 0.08 


280 


36 


77 


0.69 ± 0.27 



a Prom Gaussian fit to H I profile 
b FWHM calculated from Gaussian fit 
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